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Introduction 

It was already observed by Euler |Eull773j that the method of continued fractions 
occasionally requires a lot of tedious calculations, and even Fermat knew - as can 
be seen from the examples he chose to challenge the English mathematicians - a 
few examples with large solutions. To save work, Euler suggested a completely 
different method, which allows to compute even very large solutions of certain Pell 
equations rather easily; its drawback was that the method worked only for a specific 
class of equations. Although Euler's tricks were rediscovered on an almost regular 
basis, nobody really took this approach seriously or generalized it to arbitrary Pell 
equations. 

The main goal of |Lem2003a] and this article is to discuss certain results that 
have been obtained over the last few centuries and which will be put into a bigger 
perspective in |Lem2003b| . This is opposite to what Dickson aimed at when he 
wrote his history; in |Dicl920l vol II, preface] he says 

What is generally wanted is a full and correct statement of the 
facts, not an historians personal explanation of those facts. 

Dickson's books have occasionally been criticised for putting trivial results next 
to important ideas, that is, for not separating the wheat from the chaff; observe, 
however, that a history concentrating only on important ideas would hardly have 
mentioned most of the references on Pell's equation that are important for us. This 
is not because Dickson failed to see their importance, but because without the 
framework of a general theory they were hardly more than "mildly amusing" , as 
van der Poorten vdP2003 puts it in his review of Art2002 , the most recent paper 
containing examples of what we will call second 2-descents on certain Pell conies. 

In the follow-up |Lem2003b"| to this article, I will explain the theory of the first 
2-descent, study parts of the Selmer and Tate-Shafarevich groups attached to Pell 
conies, and interpret the results discussed here from this modern point of view. 



1. EULER 

1.1. The Content of Euler's Article. In |Eull773| . Euler remarked that the 
method for solving Pell's equation x 2 — dy 2 = 1 he has given in |Eull765] (essentially 
equivalent to the CFM) is very powerful, but that there are certain values of d 
for which even the CFM produces the solution only after tedious calculations; he 
mentions the example d = 61, where the smallest positive solution is given by 
x = 1766319049 and y = 226153980. 
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Eulcr then goes on to describe how solutions of certain auxiliary equations lead 
to solutions of the Pell equation. His first equation is 

q 2 - ap 2 = -1. (1) 

He writes 

Problema 1. Si fuerit app— 1 = qq, invenire numeros x et y, ut fiat 
axx + 1 = yy. 1 

Euler multiplies ap 2 — 1 = q 2 through by Aq 2 and adds 1 to get Aap 2 q 2 + 1 = 
Aq 4 + Aq 2 + 1 = (2q 2 + l) 2 . He has proved: 

Lemma 1.1. // — 1 = q 2 — ap 2 , then y 2 — ax 2 = 1 for x = 2pq and y = 2q 2 + 1. 

Now Euler investigates the new equation ap 2 — 1 = q 2 more closely: 

Problema 2. Investigare numeros a, pro quibues fieri potest app — 
1 = qq, hincque ipsos numeros x et y assignare, ut fiat axx+1 = yy. 2 

Euler writes ap 2 = q 2 + 1 in the form a = q t" 1 and observes that he has to find 
p and q in such a way that this fraction becomes an integer. Since a and p (and 
therefore also p 2 ) are sums of two squares, there exist integers b,c,f,g such that 
p 2 = b 2 + c 2 and q 2 + 1 = (b 2 + c 2 )(/ 2 +g 2 ), so in particular a = f 2 +g 2 . Comparing 
both sides he deduces that, for an appropriate choice of these numbers, we must 
have q — bf + eg and ±1 = bg — cf. 

Now Eulcr assigns values to b, c, f, g in such a way that ±1 = bg — cf (which can 
be done in infinitely many ways), and then computes a, p, q = bf + eg and finally 
x = 2pq and y = 2q 2 + 1. Actually, he starts by fixing p = 5, which leads to b = 3, 
c = A in view of p 2 = 25 — 3 2 + 4 2 . Then he computes the following table: 



/ 


l 


2 


4 


5 


7 


8 


9 


l 


3 


5 


7 


9 


11 


a 


2 


13 


41 


74 


130 


185 


q 


7 


18 


32 


43 


57 


68 


X 


70 


180 


320 


430 


570 


680 


y 


99 


649 


2049 


3699 


6499 


9249 



Table 1. 



He also remarks that the cases in this table are not very difficult, and gives the 
additional examples p = 13, p = 17, and p = 25. 

For ease of reference, let us collect Euler's result in the following 

Proposition 1.2. The equation —1 = q 2 — ap 2 is solvable if and only if there exist 
f,g£N with a = f 2 +g 2 such that there are b,c£N withp 2 = b 2 +c 2 , bg — cf = ±1. 
and q = bf + eg. 

Now Euler turns to the next equation, namely 

q 2 - ap 2 = -2. (2) 

He first observes f |Eull773l Problema 4]) 

^Problem 1. Given ap 2 — 1 = q 2 , to find numbers x and y such that axx + 1 = yy. 
^Problem 2. To investigate the numbers a for which we can solve ap 2 — 1 = q 2 , and then to 
assign the numbers x and y satisfying ax 2 + 1 = y 2 . 
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Lemma 1.3. If —2 = q 2 — ap 2 , then y 2 — ax 2 = 1 for x = pq and y = q 2 + 1. 

As above he then deduces from ap 2 = q 2 + 2 that a = f 2 + 2g 2 for some f,g € N, 
and that p 2 = b 2 + 2c 2 , and then concludes that cf — bg = ±1. 

Proposition 1.4. The equation —2 = q 2 — ap 2 is solvable if and only if there 
exist f, g E N with a = f 2 + 2g 2 such that there are b, c E N with p 2 = b 2 + 2c 2 , 
bg — cf = ±1, and q = bf + 2cg. 

Euler's first example is p = 3, which in view of 3 2 = l 2 + 2 • 2 2 implies 6=1 
and c = 2. Euler solves 2/ — g = ±1, giving a = f 2 + 2g 2 , q = f + Ag, hence the 
solutions x = 3q and y = q 2 + 1 of the Pell equation 1 = y 2 — ax 2 : 



f 


l 


1 


2 


2 


3 


3 


4 


9 


l 


3 


3 


5 


5 


7 


7 


a 


3 


19 


22 


54 


59 


107 


114 


q 


5 


13 


14 


22 


23 


31 


32 


X 


15 


39 


42 


66 


69 


93 


96 


y 


26 


170 


197 


485 


530 


962 


1025 



Table 2. 



In addition, Euler discusses the examples p — 9, 11, 17, 19 and then (Problema 
5) goes on to investigate the equation 

q 2 - ap 2 = +2. (3) 

Lemma 1.5. If2 = q 2 — ap 2 , then y 2 — ax 2 = 1 for x = pq and y = q 2 — 1. 

His main result in this case is 

Proposition 1.6. The equation 2 = q 2 — ap 2 is solvable if and only if there exist 
f,g E N with a = f 2 — 2g 2 such that there are b, c E N with p 2 = b 2 — 2c 2 , 
bg — cf = ±1, and q = bf — 2cg. 

As examples, Euler treats the cases p = 7, 17, 23. 

In Problema 7, Euler studies the equation app + 4 = qq: 

Lemma 1.7. 7/4 = q 2 — ap 2 , and if p and q are odd, then y 2 — ax 2 = 1 for the 

integers x = p and y = q . 

This case will not be of interest to us, so let us go right to Euler's final case, the 
equation 

s 2 - ar 2 = -4. (4) 

Lemma 1.8. If — 4 = s 2 - ar 2 , then y 2 - ax 2 = 1 for x = p and y = q 

where p = rs and q = s 2 . 

The results of Euler's problema 10 are collected in the following 

Proposition 1.9. The equation —4 = q 2 — ap 2 is solvable if and only if there 
exist f,g E N with a = f 2 — 2g 2 such that there are b, c E N with p 2 = b 2 — 2c 2 , 
bg — cf = ±1, and q = bf — 2cg. 

Euler's fourth example concerns a = 109 = 10 2 + 3 2 , 
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qui methodo vulgari molestissimos calculos requirit, 3 
Here he takes 25 2 = 625 = 24 2 + 7 2 , and finds / = 10, g = 3, b = 24 and c = 7, 
hence s = 261, giving p = 6525 and q — 261 2 + 2 = 68123 and finally 

x = 6525( 6812 f - 1 ) = 15140424455100, 

?y = 68123( 6812 f ~ 3 ) = 158070671986249. 

1.2. Interpretation with Continued Fractions. In the article discussed above, 
Euler presents solutions of the Pell equation for a variety of discriminants. The 
solutions given in TableQlcorrespond to two families, namely d = (3k— l) 2 + (4fc— l) 2 
and d = (3k + l) 2 + (4k + l) 2 . Developing \fd into continued fractions (which Euler 
did not do) we find: 



d 


Vd 


(3k- 


I) 2 


+ (4k 


-l) 2 


= 5 2 fc 2 


- 14fc + 2 


[5fe - 2, 1, 1, 1, 1, lQfc- 4j 


(3k + 


I) 2 


+ (4k 


f l) 2 


= 25k 2 


+ 14fc + 2 


[5k + 1, 2, 2, lOfc + 2] 


(5k -2] 


2 + 


(12k- 


5) 2 = 


-- 13 2 k 2 


- UOk + 29 


[13fc-6,l,l,l,l,l,l,26fc- 12] 



The last line comes from Euler's second table, which we did not reproduce here. 

In general, assume that (r, s, t) is a primitive Pythagorean triple with s even. 
Then there exist m, n such that r = m 2 — n 2 , s = 2mn and t = m 2 + n 2 . Euler has 
to solve the linear equation (m 2 — n 2 )f — 2mng = ±1. 

Let us look at families that admit n = 1 as a solution. Then we have to consider 
(m 2 - 1)/ - 2mg = ±1; if we put m = 2k, we find (4k 2 - 1)/ - 4kg = ±1. The 
solutions of this equation are / = 1 + 4ku, g = k + (4k 2 — l)u and / = — 1 + 4ku, 
g = —k + (4k 2 — l)u for u = 0, 1, 2, As above we find 



d 


Vd 


[(4fc 2 + l)w + fc] 2 +4fcu+l 


[e, 2k, 2k, 2e] 


e = (4fc 2 + l)u + k 


[(4fc 2 + l)w-fc] 2 -4fcu+l 


[e, 1,2/5+ 1,2k + 1,1, 2e] 


e = (4k 2 + l)u-k-l 



2. Hart 

2.1. Hart's Article. In a short note Harl878a , Hart studied the solvability of 
the diophantine equation X 2 — dY 2 = —1; let us quote the first half of Harl878a : 
To find general values of x and y to solve the problem 

x 2 - Ay 2 = -1. 

Let A be a non-quadrate number = the sum of two squares = r 2 +s 2 , 
then we shall have 

2 / 2 i 2\ 2 -i 

x — (r + s )y = —1, 

and by transposition 

x 2 - r 2 y 2 = s 2 y 2 - 1 = (sy - l)(sy + 1); 

whence 

x 2 = r 2 y 2 + (sy-l)(sy + l)=D; .\ x = [r 2 y 2 + (sy - l)(sy + 1)] 1/2 . 

''for which the usual method requires the most tedious calculations. 
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Let r 2 y 2 + (sy — l)(sy + 1) = [ry — (m 4- n)(sy — l)] 2 . Reducing 
this we get 

m 2 + n 

y = — 



■ — 2nrm — sn 2 

Here, in order to have y integral, put sm 2 — 2nrm — sn 2 = ±1; 
whence, transposing and dividing by s, we have 

2nr sn 2 ± 1 



m 2 m 



s 

and by quadratics 



rn ± \l (r 2 + s 2 )n 2 ± s „ 

m= — , .'. y = ±(m 2 +n 2 ). 

s 

In the general value of m, r and s may represent any numbers one 
of which is even, and the other any odd number except 1, and n 
can be found by trial, or, if large, by the solution of the formula 

Let us now give a modern interpretation of Hart's idea. Assume we want to solve 

x 2 - Ay 2 = -1 (5) 
for some squarefree integer A. It was already known to Brahmagupta (see Whitford 
Whil912 J that this implies that A is the sum of two squares. Thus A = r 2 + s 2 for 
integers r, s, and we have to solve ir 2 = r 2 y 2 + (sy — l)(sy + l). Now parametrize this 
conic using the rational point P — (x, y) — (J, -): the equation x — ry — t(sy — 1) 
describes lines through P with slope ;qr^-; if we pick t = ^ rational, such a line 
will intersect the conic in another rational point, and a simple calculation shows 

2 , 2 

that its w-coordinate is y — — $ + " r . We want values m and n for which y is 

integral; thus we are led to consider 

sm 2 — 2nrm — sn 2 = ±1. (6) 

Any integral solution of 10 will give an integral solution of JSJ). Euler and Lagrange 
have shown how to reduce © to a Pell equation: interpreting © as a quadratic 
polynomial in m, a necessary condition for the existence of an integral solution is 
that the discriminant 4r 2 n 2 + 4s(sn 2 =F 1) be a square, i.e., that one of the equations 

P 2 - An 2 = ±s (7) 
have an integral solution. This shows 

Proposition 2.1. Let d = r 2 + s be a sum of two squares; if the equation sm 2 — 
2rmn — sn 2 = ±1 has an integral solution (m,n), then so does x 2 — dy 2 = —1; in 
fact, we can put y = m 2 + n 2 . 

Thus Hart has proved that the solvability of © for some choice of r, s implies 
the solvability of it does not follow (at least not directly) from Hart's proof that 
the condition is also necessary. The necessity actually was proved by Euler: Hart's 
result is nothing but a special case of Propositior ll.21 In fact, if JSJ) is solvable, then 
by Euler there exist r, s with A = r 2 + s 2 such that y 2 — b 2 + c 2 and bs — cr = ±1. 
Since (b,c,y) is a Pythagorean triple, we can write b = m 2 — n 2 , c — 2mn and 
y = m 2 + n 2 ; plugging the values of b and c into bs — cr — ±1 then yields 10- Thus 
we have shown that combining the results by Euler and Hart gives 
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Proposition 2.2. The equation x 2 — Ay 2 = — 1 is solvable if and only if there exist 
r, s € N with A = r 2 + s 2 such that the diophantine equation Jj^) has an integral 
solution. 

2.2. The Negative Pell Equation. Hart's equations were mentioned by Dickson 
|Dicl920j and rediscovered by Sansone |Sanl925allSanl925b| and Epstein |Epsl934| ; 
these two authors used the fact that solvability of the negative Pell equation implies 
the solvability of x 2 — dy 2 = a (Hart's equation Q), where d = a 2 + b 2 and a is 
odd. Reducing this equation modulo primes p dividing d shows that (a/p) — +1: 

Proposition 2.3. Let d be a squarefree integer. If x 2 — dy 2 = —1 is solvable, then 
d = a 2 + b 2 with a odd, and (a/p) = +1 for all primes p | d. 

Escott Escl905 had already proved that if x 2 — Dy 2 = — 1 is solvable in integers 
and D = a 2 + b 2 , then (a/D) = 1 or (b/D) — 1. He also showed that the condition 
is not sufficient by using the example 2306 = 41 2 + 25 2 . Recently, Proposition ^. 31 
was rediscovered by Khessami Pilerud KP1999 . 

Grytczuk, Luca, & Wojtowicz GLW2000 proved the following result, which is 
just a slightly reformulated version of Euler's Proposition ^21 

Proposition 2.4. The equation s 2 — dr 2 = — 1 has an integral solution if and only 
if there is some odd integer B G Z with d = A 2 + B 2 and a Pythagorean triple 
(a, 6, c) such that a A — bB = ±1. 

3. Sylvester 

Sylvester |Syll88l| proved the following 

Proposition 3.1. Let A = 2f 2 + g 2 be a prime with f odd; then the equation 
fy 2 + 2gxy — 2 fx 2 = ±1 is solvable in integers. 

For the proof, let (u, v) be the minimal solution of the Pell equation u 2 ~ Av 2 = 1. 
Playing the usual game, Sylvester arrives at p 2 — Aq 2 — 1 or p 2 — Aq 2 — —2, the 
other signs being excluded because of A = 3 mod 8. The minimality of the solution 
implies that we have p 2 — Aq 2 = — 2, and using unique factorization in Z[ v / — 2 ] we 
get p + v 7 — 2 = (g + f v 7 — 2 )(y + x\J — 2 ) 2 . Comparing the imaginary parts shows 
that fy 2 + 2gxy - 2 fx 2 = ±1. 

Sylvester's result can be derived from Euler's Proposition !! ,4l in the same way we 
deduced Hart's result from Proposition ll.2l all we have to do is replace Pythagorean 
triples by solutions of the equation x 2 + 2y 2 = z 2 . 

4. Gunther 

In |Guel882 * | . S. Gunther discussed the comments of Theon Smyrnaeus on Plato's 
work and concluded that he must have been familiar with the Pell equation 

2x 2 -l = y 2 . (8) 

He then tries to reconstruct a possible approach to the solution of this equation. 

Gunther suggests writing (JSJ) as x 2 — 1 = y 2 — x 2 , substitutes x + 1 = |(y + x) 
and x — 1 = ^ (y — x) , then sets 

p 2 + 2pq -q 2 = z, (9) 
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solves this equation for p, and then deduces 

4q 2 + z =F 2qJ2q 2 + z -4g 2 -z± ^qj2q 2 + z 
x= , y = . 

z z 

Thus any solution of © with z = ±1 will lead to an integral solution of the negative 
Pell equation (JSJ . 

The same method, he remarks, works for the more general equation 

(a 2 + b 2 )x 2 -l = y 2 , 

which can be written in the form (ax — l)(ax + 1) = (y — bx)(y + bx). Substituting 
ax+l = £(y + bx) and ax—1 = -(y—bx) leads to equations equivalent to Hart's © 
and J7J), but the formulas derived by Giinther are incorrect, and his final conclusion 
is unclear. 

5. Gerardin 

While extending existing tables of solutions of the Pell equation by Legendre, 
Bickmore, and Whitford, A. Gerardin |Gerl917| complained about the tedious work 
necessary when using the theory of continued fractions: 

La recherche pratique de la solution minima etait faite jusqu'a 
present sur les fractions continues, ce qui demande en general beau- 
coup de soins et de temps. 

He then presents, by giving a few examples, a 'new method' for solving Pell equa- 
tions; he uses the auxiliary equations 

x 2 -Ay 2 = ±4, ±2,-1 (10) 

and remarks that from their solutions one can pass easily to the solution of the 
corresponding Pell equation. 

For A = 941 = 29 2 + 10 2 he solves the pair of equations 

31 2 - 941 • l 2 = +2- 10, 
184 2 - 941 • 6 2 = -2 • 10, 

and then constructs a solution of x 2 — 941y 2 ~ —4 by setting y = l 2 + 6 2 = 37, 
giving x = 1135. 

The general cases are given by the following formulas: 

a) z 2 - At 2 = -1, A = m 2 +n 2 ,t = a 2 + (3 2 

(ma - nf3) 2 - Af3 2 = ±m, (11) 
(mj3 + no) 2 - Aa 2 = ^m. (12) 

b) z 2 - At 2 = +2, A = m 2 - 2n 2 , t = a 2 - 2/3 2 

(not - m/3) 2 - A0 2 = ±n, (13) 
(n(3 - ma) 2 - Aa 2 = ±2n. (14) 



The practical search for the minimal solution was made up until now via continued fractions, 
which demands in general a lot of care and time. 
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c) z 2 - At 2 = -4, A = m 2 +n 2 , t = a 2 + (3 2 

(na - m/3) 2 - A(3 2 = ±2n, (15) 
(n/3 + ma) 2 - Aa 2 = T^n. (16) 

Gerardin also remarks that the case z 2 — At 2 — — 2 can be treated similarly. 

Thus in case a), Gerardin writes A = m 2 + n 2 and then tries to solve the two 
equations r 2 — As 2 — m and t 2 — Au 2 = —m; putting t — a 2 + (3 2 then gives a 
solution of z 2 — At 2 = — 1, from which a solution to the Pell equation X 2 — AY 2 = 1 
is easily derived. Observe the similarity with the result of Hart |Harl878a] discussed 
above. 

Gerardin does not give any proofs, but his claims are easily verified. Let us first 
consider equation jH) . Plugging in A = m 2 + n 2 and simplifying we get 

ma 2 - 2na/3 - mf3 2 = ±1. (17) 

Observe that taking \12\ would lead to the very same equation; in particular, i|ll|) 
and (|12l) are equivalent. Also note that m must be odd for (|17fl to be solvable. A 
simple calculation now shows that 

{na 2 - 2ma(3 - n(3 2 ) 2 - A(a 2 + 1 ) 2 = 

-(ma 2 - 2na/3- m/3 2 ) 2 = -1. (18) 

Similarly, in cases b) and c) we get the equations 

no 2 - 2maf3 + 2n/3 2 = ±1 (19) 
na 2 - 2maf3 - n/3 2 = ±2 (20) 

(it is easy to see that, in H20JI . m must be odd) as well as 

(ma 2 - Ana/3 + 2m0 2 ) 2 - A(a 2 - 2/3 2 ) 2 = 

2(na 2 - 2maf3 + 2nf3 2 ) 2 =+2 (21) 
(ma 2 - 2na(3 + m(3 2 ) 2 - A(a 2 + f3 2 ) 2 = 

-(no 2 - 2maj3 - n(3 2 ) 2 =-4. (22) 

Comparing (|10|) with Euler's equations CEJ - Q from Section ^ it will come as 
no surprise that Gerardin's formulas are an easy consequence of Euler's results. 

6. Hardy & Williams, Bapoungue, Arteha 

6.1. Hardy & Williams. In |HW1986| . K. Hardy & K. Williams investigate the 
solvability of the diophantine equation 

fx 2 -2gxy~fy 2 = 1. (23) 

Their main result is 

Proposition 6.1. The equation —1 = q 2 — ap 2 is solvable if and only if there exist 
/, g € N with a = f 2 + g 2 such that there are x, y £ Z with \2!$ . In this case, the 
pair (/, g) with f odd is unique. 
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Apart from the uniqueness assertion, this result is an almost trivial conse- 
quence of Euler's Proposition 11.21 we start with the observation that (b,c,p) is 
a Pythagorean triple. Now — 1 = q 2 — ap 2 implies that a and p are odd. Assuming 
that c is odd, we find from bg — cf — ±1 that / is also odd. Changing the signs 
of /, g if necessary we may assume that bg — cf = — 1 . The parametrization of 
Pythagorean triples shows that b = 2xy and c = x 2 — y 2 ; the equation bg — cf = —1 
then becomes f(x 2 — y 2 ) — 2gxy = 1. This proves Proposition 16. II except for the 
uniqueness part. 

The claim that there is essentially only one such pair (/, g) is an important 
contribution: as we will see later, it should be seen as (part of) an analogue of 
Dirichlet's Theorem |Lem2003al Thm. 3.3]. 

6.2. Bapoungue. Inspired by the work of Hardy & K. Williams HW1986 , Ba- 
poungue |Bapl989[ |Bapl998| |B"ap2000a| |Bap"2000b| |Bap2002| started investigating 
the solvability of the diophantine equation 

ax 2 + 2bxy - kay 2 = ±1 (24) 

for values of k for which k ) has class number 1. The identity 

(-bx 2 + 2akxy + kby 2 ) 2 - {b 2 + ka 2 )(x 2 + ky 2 ) 2 = -k(ax 2 + 2bxy - kay 2 ) 2 

yields, upon substituting a solution of (|24Jl for (x 7 y), the equation 

{-bx 2 + 2akxy + kby 2 ) 2 - S(x 2 + ky 2 ) 2 = -k. (25) 

This shows ( |Bapl998| Thm. 2]) 

Theorem 6.2. If \2$ has an integral solution, then so does 

Multiplying l|24|l through by a and completing the square we get 

(ax + by) 2 - 5y 2 = ±a, where 5 = ka 2 + b 2 . (26) 

Thus the solvability of 1)24(1 implies the solvability of the Pell equation X 2 — SY 2 = a, 
where solvability denotes solvability in integers. 

Similarly, multiplying l)24|l through by —ka and completing the square we get 

{kay - bx) 2 - Sx 2 = Tka. (27) 

The special cases k = 1, 2 of these equations go back to Euler, Hart, and Gerardin. 
The main result of Bapoungue's thesis |Bapl989~| is 

Theorem 6.3. Let k <E {2,3,7,11,19,43,67,163} (this implies that Q(y/^k) has 
class number 1 ). If a, b are positive integer with a odd such that p = ka 2 + b 2 is 
prime, then \2$ is solvable if and only if 12 is solvable. 

A similar result holds for k = 1. The case k = 2 is Sylvester's Proposition ^. II 
In |Bap2000b| , the following result is proved: 

Theorem 6.4. Let k be as above. Among all pairs (a, b) of coprime natural numbers 
with a odd and d = ka 2 + b 2 , there is exactly one pair for which \2J$ is solvable. 

This generalizes Theorem 16. II 
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6.3. Arteha. The last rediscovery of the method of Euler-Hart-Gerardin so far is 
due to Arteha Art2002 ; his results are 

Proposition 6.5. Consider the Pell equation 

x 2 -dy 2 = 1 (28) 

for primes d. 

(1) If d = 1 mod 4, and write d = a 2 +b 2 with a odd. Then the minimal positive 
solution of Pell's equation $28\) is given by 

y = 2\2arrm + b(m 2 — n 2 )\(m 2 + n 2 ), 

where m and n satisfy 

a(m 2 — n 2 ) — 2bmn = ±1. 

(2) If d = 3 mod 8, write d = a 2 + 2b 2 . Then the minimal positive solution of 
Pell's equation \2fy) is given by 

y = \ibmri + a\m 2 — 2n 2 \ | (m + 2n 2 ), 

where m and n satisfy 

b\m 2 - 2n 2 \ - 2amn = ±1. 

(3) If d = 7 mod 8, write d = a 2 — 2b 2 . Then the minimal positive solution of 
Pell's equation is given by 

y = \(a(m 2 + 2n 2 ) - 4bmn)(m 2 - 2n 2 )\, 

where m and n satisfy 

2amn — b(m 2 + 2n 2 ) = ±1. 

7. Summary 

Starting with Euler in |Eull773j . many authors have come up with essentially 
the same idea: solving the Pell equation 

X 2 - dY 2 = 1 (29) 

becomes easier by looking at certain auxiliary equations. 
The first step is writing down Legendre's equations 

rx 2 - sy 2 = 1,2 for d = rs. (30) 

There is one nontrivial equation among these with a solution, and the smallest 
solution will have about half as many digits as the smallest solution of (|29l) . 

The second step is to look at equations whose solutions give rise to a solution 
of one of the equations l|3"0"|) : but this second step has never been completed in 
full generality. What we have are specific equations applicable only in special 
situations; these were first discovered by Euler, rediscovered by Hart, Sylvester, 
Gunther, Gerardin, Hardy & Williams, and Arteha, and slightly generalized by 
Bapoungue. 



HIGHER DESCENT ON PELL CONICS 



11 



References 

[Art2002] S.N. Artoha, Method of hidden parameters and Pell's equation, JPJ Algebra Number 

Theory Appl. 2 (2002), 21-46; cf. p. ITllTol 
[Bapl989] L. Bapoungue, Sur la resolubilite de I'equation ax 2 + 2bxy — kay 2 = ±1, These Univ. 

Caen, 1989; see also C. R. Acad. Sci. Paris 309 (1989), 235-238; cf. p. El 
[Bapl998] L. Bapoungue, Un critere de resolution pour I'equation diophantienne ax 2 + 2bxy — 

kay 2 = ±1, Expos. Math. 16 (1998), 249-262; cf. p. El 
[Bap2000a] L. Bapoungue, Sur la resolubilite de I'equation ax 2 + 2bxy — Say 2 = ±1, IMHOTEP, 

J. Afr. Math. Pures Appl. 3 (2000), 97-111; cf. p. El 
[Bap2000b] L. Bapoungue, Sur les solutions generates de I'equation diophantienne ax 2 + 2bxy — 

kay 2 = ±1, Expos. Math. 18 (2000), 165-175; cf. p. El 
[Bap2002] L. Bapoungue, The diophantine equation ax 2 + 2bxy — Aay 2 = ±1, Intern. J. Math. 

Math. Sci. 35 (2003), 2241-2253 El 
[Dicl920] L.E. Dickson, History of the Theory of Numbers, vol I (1920); vol II (1920); vol III 

(1923); Chelsea reprint 1952; cf. p. Q]EI 
[Epsl934] P. Epstein, Zur Aufldsbarkeit der Gleichung x 2 — Dy 2 = —1, J. Reine Angew. Math. 

171 (1934), 243-252; cf. p. El 
[Escl905] E.B. Escott, Solution de I'equation x 2 - Dy 2 = -1, LTntermed Math. 12 (1905), 53; 

cf. p. El 

[Eull765] L. Euler, De usu novi algorithmi in problemate Pelliano solvendo, Novi Acad. Sci. 

Pctropol. 11 (1765) 1767, 28-66; Opera Omnia 1-3, 73-111; cf. p. □ 
[Eull773] L. Euler, Nova subsidia pro resolutione formulae axx + 1 = yy, Sept. 23, 1773; Opusc. 

anal. 1 (1783), 310; Comm. Arith. Coll. II, 35-43; Opera Omnia 1-4, 91-104; cf. p. fTH^llTTTl 
[Gerl917] A. Gerardin, Sur I'equation x 2 - Ay 2 = 1, L'Ens. math. 19 (1917), 316-318; Sphinx- 

(Edipe 12 June 15, 1917, 1-3; cf. p. □ 
[GLW2000] A. Grytczuk, F. Luca, M. Wojtowicz, The negative Pell equation and Pythagorean 

triples, Proc. Japan Acad. 76 (2000), 91-94; cf. p. El 
[Guel882] S. Giinther, Ueber einen Specialfall der Pell'schen Gleichung, Blatter fur das Bay- 

erische Gymnasial- und Realschulwesen 17 (1882), 19—24; cf. p. El 
[HW1986] K. Hardy, K. Williams, On the solvability of the diophantine equation dV 2 — 2eVW — 

dW 2 = 1, Pac. J. Math. 124 (1986), 145-158; cf. p. E1E1 
[Harl878a] D.S. Hart, Solution of an indeterminate problem, Analyst 5 (1878), 118-119; cf. p. El 

[KP1999] Kh. Khessami Pilerud, On the Diophantine equation x 2 — Ny 2 = — 1, (Russian) Vestnik 

Moskov. Univ. Ser. I Mat. Mekh. (1999), no. 2, 65-67; Engl, transl. Moscow Univ. Math. Bull. 

54 (1999), no. 2, 48-49; cf. p. El 
[Lem2003a] F. Lemmermeyer, Higher Descent on Pell Conies I. From Legendre to Selmer, 

preprint 2003; cf. p. E1E1 
[Lem2003b] F. Lemmermeyer, Higher Descent on Pell Conies III. The First 2-Descent, preprint 

2003; cf. p. □ 

[vdP2003] A. van der Poorten, Review 2003i:11040, MathSciNet; cf. p. Q] 

[Sanl925a] G. Sansone, Sulle equazioni indeterminate delle unitd di norma negativa dei corpi 
quadratici reali, Rend. Acad. d. L. Roma (6) 2 (1925), 479-484; cf. p. El 

[Sanl925b] G. Sansone, Ancora sulle equazioni indeterminate delle unitd di norma negativa dei 
corpi quadratici reali, Rend. Acad. d. L. Roma (6) 2 (1925), 548-554; cf. p. El 

[Syll881] J.J. Sylvester, Mathematical Question 6243, Educational Times 34 (1881), 21-22; cf. 

p. El 

[Whil912] E.E. Whitford, The Pell equation, New York 1912, 193 pp; cf. p. El 



Department of Mathematics, Bilkent University, 06800 Bilkent, Ankara, Turkey 
E-mail address: franzafen.bilkent.edu.tr 



